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Current cosmological constraints from Cosmic Microwave Background (CMB) anisotropics are 
typically derived assuming a standard recombination scheme, however additional resonance and ion- 
izing radiation sources can delay recombination, altering the cosmic ionization history and the cosmo- 
logical inferences drawn from CMB data. We show that for recent observations of CMB anisotropy, 
from the Wilkinson Microwave Anisotropy Probe satellite mission 5-year survey (WMAP5) and from 
the ACBAR experiment, additional resonance radiation is nearly degenerate with variations in the 
spectral index, tIs, and has a marked effect on uncertainties in constraints on the Hubble constant, 
age of the universe, curvature and the upper bound on the neutrino mass. When a modified recom- 
bination scheme is considered, the redshift of recombination is constrained to z, = 1078 ± 11, with 
uncertainties in the measurement weaker by one order of magnitude than those obtained under the 
assumption of standard recombination while constraints on the shift parameter are shifted by Icr 
to 72. = 1.734 ± 0.028. Although delayed recombination limits the precision of parameter estima- 
tion from the WMAP satellite, we demonstrate that this should not be the case for future, smaller 
angular scales measurements, such as those by the Planck satellite mission. 
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I. INTRODUCTION 

The recent measurements of the Cosmic Microwave 
Background (CMB) flux provided by the five year 
Wilkinson Microwave Anisotropy Probe (WMAP) 
mission (see [l|, [2| and the ACBAR collaboration (see 
Q) have confirmed several aspects of the cosmolog- 
ical standard model and improved the constraints 
on several key parameters. The constraints on the 
neutrino mass, for example, coming just from CMB 
data are drastically improved and are now competi- 
tive with current direct laboratory measurements (see 
e-g- 0, [j])- Moreover, the presence of a neutrino 
background is now inferred at more than 95% c.l. by 
CMB data alone (see e.g., @,[1|)- Finally, infiationary 
parameters, curvature, baryonic and dark matter den- 
sities are also now better determined due to improved 
treatment of systematics (see 0, @)- 

These spectacular results, apart from the experi- 
mental improvements, have been possible due to the 
high precision of the CMB theoretical predictions that 
have now reached an accuracy close to 0.1% over a 
wide range of scales. A key ingredient in the CMB 
precision cosmology is the accurate computation of 
the recombination process. Since the seminal papers 
by Peebles and Z'eldovich (see [3, Q) detailing the re- 
combination process, further refinements to the stan- 
dard scheme were developed [91] allowing predictions 
at the accuracy level found in data from the WMAP 
satellite and predicted for the future Planck satellite 

While the attained accuracy on the recombina- 
tion process is impressive, it should be noticed that 
these computations rely on the assumption of stan- 
dard physics. Non-standard mechanisms such as (just 
to name a few) high redshift stars or active galac- 
tic nuclei, topological defects and dark matter decays 
could produce extra sources of radiation and modify 



the recombination process. With the WMAP results 
and the future Planck data, it therefore becomes con- 
ceivable that deviations from standard recombination 
may be be detected. 

Several papers in recent years have indeed inves- 
tigated this possibility. For example, one could use 
a phenomenological approach such as in [lj| or |14| . 
Other works have adopted a more physically moti- 
vated (but model-dependent) approach and consid- 
ered modified recombination by allowing time vari- 
ations in the fine-structure ([l5|) or gravitational con- 
stants (lial)- Here we instead focus on delayed re- 
combination mechanisms based on the hypothesis of 
extra sources of ionizin g a nd resonance radiation at 
recombination (see e.g. [ITI [la. [l9| ) . 

Extra photon sources can be generated by a vari- 
ety of mechanisms. A widely considered process is 
the decay or annihilation of massive particles [20, [2lj, 
m, [H, [ll, [H, [11]. The decay channel depends on 
the nature of the particles, and could, for example, 
include charged and neutral leptons, quarks or gauge 
bosons. These particles may then decay further, lead- 
ing to a shower cascade that could, among other prod- 
ucts, generate a bath of lower energy photons that 
would interact with the primordial gas and cosmic 
microwave background. Interestingly these models, 
as well as injecting energy at recombination, boost 
the ionization fraction after recombination and can 
distort the ionization history of the universe at even 
later times, du ring galaxy formation and reionization 
[27, 23^, 2^ [30; [m ■ Other mechanisms include evap- 
oration of black holes [1^, [84] or inhomogenities in 
baryonic matter [l8j . 

Several authors have already compared delayed re- 
combination with cosmological data (see for example 
m, [H, [H, [H, [13]). With respect to these previ- 
ous analyses, we assess below the improvements given 
by more recent data from the WMAP five year sur- 



vey and from the ACBAR experiment. Moreover, we 
will study in detail the impact of delayed recombi- 
nation on the current constraints on cosmic param- 
eters. As showed in [l7| delayed recombination has 
two main effects: damping of the CMB anisotropy 
and polarization at small angular scales and a shift of 
the acoustic peaks in their angular spectra. As wc will 
see, these two effects change in a significant way the 
current constraints on such parameters as the spectral 
index ng and the Hubble constant Hq. We then study 
the impact on the current determination of the red- 
shift of recombination z* and on the shift parameter 
7?.. Both parameters are used to provide complemen- 
tary geometric constraints on the background expan- 
sion history; i? as a measure of the angular diameter 
distance to last scattering, and z* in interpreting the 
scale of baryon acoustic oscillations \38, 6^ , which in 
combination with supernova type la data, can impose 
constraints on the dark energy component (see e.g. 
[40l . l4ll . |42| . Here we show that delayed recombina- 
tion has decisive effects on the determination of these 
parameters. Finally, we discuss the implications of a 
modified recombination scheme for the Planck satel- 
lite. 

The paper proceeds as follows: in section |TT] we 
briefly describe the delayed recombination scheme. In 
lIHI we analyze the latest CMB data and place new con- 
straints on delayed recombination. Moreover, we fore- 
cast the constraints prospectively attainable by the 
Planck satellite mission, and study the impact that 
a modified recombination scheme can have on several 
key cosmological and astrophysical parameters. In llVI 
we draw our conclusions and review their future im- 
plications. 



II. A MODIFIED IONIZATION HISTORY 

Following [3, 3 we can model the evolution of the 
electron ionization fraction, Xe in a simplified manner 
for the recombination of hydrogen: 
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where n is the number density of atoms, Oc and be 
are the effective recombination and photo-ionization 
rates for principle quantum numbers > 2, AB is the 
difference in binding energy between the 1*** and 2"'' 
energy levels and 
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where Aq is the wavelength of the single Ly-a transi- 
tion from the 2p level, Ais2s is the decay rate of the 
metastable 2s level, nis = n(l — Xe) is the number of 
neutral ground state H atoms, and H{z) is the Hubble 
expansion factor at a redshift z. 



As in [3J|, we include the possibility of extra pho- 
tons at key wavelengths that could modify this recom- 
bination picture, namely resonance (Ly-a) photons 
with number density, tIq, which promote electrons to 
the 2p level, and ionizing photons, n, . \VA UM. [l9l [SSf 

—JT = £a{z)H{z)n, —^=£,(z)H{z)n. (3) 
at at 

This leads to a modified evolution of the ionization 
fraction 

We employ the widely used RECFAST code , in the 
cosmomc package [43 modifying the code as in ([4]) to 
include two extra constant parameters, £„ and e.j. 

The effects of delayed recombination on the primor- 
dial power spectra have been througly investigated by 
jl7| . Namely, the main consequences are damping 
of CMB anisotropics and polarization on small an- 
gular scales and a shift of the acoustic peaks towards 
larger scales. This introduces a main degeneracy with 
two parameters: the scalar spectral index ng and, in 
a flat universe, the Hubble constant Hq. In Figure 
[T] we indeed plot anisotropy, polarization and cross 
temperaturc-anisotropy power spectra for two degen- 
erate models. As we see, when delayed recombination 
is included, the degenerate spectra have a larger spec- 
tal index to compensate the damping and a smaller 
Hubble parameter to compensate the shift. This de- 
generacy will clearly appear as a correlation between 
the constraints on these parameters as we show in the 
next section. 

We compare delayed recombination with current 
data by making use of the the publicly available 
Markov Chain Monte Carlo package cosmomc |43j . 
Other than Eq. and e^, we sample the following seven- 
dimensional set of cosmological parameters, adopting 
flat priors on them: the physical baryon and CDM 
densities, uji, ~ Vliih? and Wc = fig^^, the Hubble pa- 
rameter iJo, the scalar spectral index, rig, the nor- 
malization, h\\{)^^As{k = 0.05/Mpc), and the optical 
depth to rcionization, r. We also consider the pos- 
sibility of having a curved universe with Uk ^ [44| 
and three massive neutrinos each with same mass vfii, 
and overall energy density: 
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Wc consider purely adiabatic initial conditions. In 
what follows we also consider bounds on the derived 
"shift" parameter TZ defined as (see j40|): 
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FIG. 1: Best fit WMAP5 Spectra and degenerate 
Anisotropy (Top), Polarization (Middle Panel) and cross 
Anisotropy-Polarization (Bottom Panel) angular power 
spectra in case of delayed recombination. 



where sinn(a::) = {sin(a;),a:,sinh(a;)} for closed, flat 
and open geometries respectively, with 
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a, V^mfl + t^kfl^ + ti^Aa^ + UJ\ 
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where uJm = ujb + ^c and a is the scale factor with 
a* = (1 + z*)~^. The MCMC convergence diagnos- 
tic tests are performed on 8 chains using the Gclnian 
and Rubin "variance of chain mean" / "mean of chain 
variances" R statistic for each parameter. Our 1 — D 
and 2 — D constraints arc obtained after marginahza- 
tion over the remaining "nuisance" parameters, again 
using the programs included in the cosmomc package. 
We use a cosmic age top-hat prior as 10 Gyr < to ^ 20 
Gyr. We include the five-year WMAP data [I (tem- 
perature and polarization) with the routine for com- 
puting the likelihood supplied by the WMAP team 
(we will refer to this analysis as WMAP5). Together 
with the WMAP data we also consider the small-scale 
CMB measurements of ACBAR Q (we will refer to 




FIG. 2: Likelihood contours at 68% and 95% c.l. in 
the 71s — Ea, plane for WMAP5 (empty contours) and 
WMAP5-I-ACBAR (filled contours) experiments respec- 
tively, allowing both ionizing and resonance radiation 
modifications to the recombination scheme. 



this analysis as WMAP5-h ACBAR). 



III. LIKELIHOOD ANALYSIS AND RESULTS 

A. Modified Recombination w^ith standard 
A-CDM 

We first analyze the effects of modified recombina- 
tion on the cosmological constraints obtained under 
the assumption of a flat, A-CDM model with mass- 
less neutrinos. After marginalization over the nui- 
sance parameters described in the previous section, 
the WMAP5 data gives £„ < 0.39 and e,, < 0.058 at 
95% c.l. {e^ < 0.31 and e^ < 0.053 when the ACBAR 
data is also included). 

In Figure [H and Figure [3] we plot the 68% and 
95% c.l. likefihood contours on the n^ — £„, and 
n, -e, planes for the WMAP5 and WMAP5-h ACBAR 
datasets respectively. The figures highlight that even 
with the improved small scale CMB data, modified 
recombination's main effect is to drive Us to higher 
values, reconciling the data with a n^ — 1, Harrison- 
Z'eldovich (HZ), spectrum, as was seen with the 
WMAP 3-year data [3J| . Any conclusion on the com- 
patibility of a particular infiationary model with the 
WMAP5 data, therefore, is highly dependent on the 
assumptions made on the recombination process and 
should be discussed with some caution in light of this. 

Marginalizing over the recombination parameters 
we indeed get n, = 0.993J:oolo (WMAP5 alone) and 
n, = 0.996+""'^2 
i.e 



- 993+0045 

- U.»»0„Q Q4Q 

_o(,38 (WMAP5-f ACBAR) at 95% c.l.., 
with Tig = 1 perfectly compatible with the data. 



Those results should be compared with the constraints 
n, = 0.959to.o27 and n, = 0.965to:o28 (95% c.l.) ob- 
tained using the same datasets and priors but with 
standard recombination. 

As we can see from Figure [2] it is possible to have 
Tis = 1 even if Cq, = 0. However this agreement, ap- 
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FIG. 3: Likelihood contours at 68% and 95% c.l. in 
tlie Us — ei, plane for WMAP5 (empty contours) and 
WMAP5+ACBAR (filled contours) experiments respec- 
tively, allowing both ionizing and resonance radiation 
modifications to the recombination scheme. 




FIG. 4: Likelihood contours at 68% and 95% c.l. on 
the Ua — ea, plane for WMAP5 assuming variations in ei 
(empty contours) and fixing e^ = (filled contours). 



parcntly in contrast with the standard WMAP5 re- 
sult, is due to the marginalization over ti which is 
still present as a free parameter in the analysis. Since 
it may be possible to have only resonance radiation at 
recombination we analyze the effects of just varying 
ta while keeping e^ = 0. In Figure [4] we plot the con- 
straints on the Us-ta plane with and without variation 
in ti. As we can see, when e^ = 0, the degeneracy be- 
tween Us and €a is more evident and HZ spectra are in 
agreement at 1-ct level with WMAP5 only if delayed 
recombination is present. 

The degeneracy with the spectral index can be eas- 
ily understood: delayed recombination suppresses the 
amplitude of the acoustic peaks in the CMB angular 
temperature anisotropy in a very similar way to an 
increase in the optical depth [l3|- This effect can be 
compensated by an increase in ns- Ionizing rcioniza- 
tion also produces significant large angular scale po- 
larization, and WMAP EE and TE constraints help 
to break the Us — ti degeneracy. However the same is 
not true for resonance recombination |3j| , and so the 
ris — ta degeneracy remains present despite improved 
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FIG. 5: Likelihood contours at 68% and 95% c.l. on the 
Ho — Eq , (top panel) and to — fa plane for WMAP5 fixing 
ei = (bottom panel). The Hubble constant is in units of 
km/s/Mpc. 



polarization measurements. 

Another effect of delayed recombination on the 
CMB spectrum, through the presence of resonance 
radiation, is a shift of the acoustic peak spectrum to- 
wards larger angular scales. An increase in e^ shifts 
the epoch of recombination towards smaller z* . This 
increases the size of the acoustic horizon at recombi- 
nation and therefore shifts the peak towards smaller 
ts (see [131: 123 )• I'^ ^ fl^t universe, this effect can 
be counter-balanced by a decrease of the Hubble con- 
stant. In Figure [5] (Top Panel) we plot the likelihood 
contours from WMAP5 on the Ho—ea plane. A degen- 
eracy with the Hubble parameter is evident and the 
analysis yields Ho = GS.llg'g at 95% c.l. that should 
compared with the value Hq ~ Tl.Oits^ at 95% c.l. as- 
suming standard recombination. Since smaller values 
of the Hubble constant are in agreement with larger 
€a it is easy to predict a correlation with the current 
age of the universe io • 

In the bottom panel of Figure [5] we plot the likeli- 
hood contours in the to — e^ plane. As we see, larger 
values of Eq are indeed in agreement with larger to- 
We find from WMAP5 that <o = 14.0+^j Gyrs to 
be compared with the constraint to = IS.Til^o'g ^^~ 
der standard recombination. Current CMB age con- 
straints can, therefore, be underestimated if one as- 
sumes a standard recombination scheme. 

The impact of delayed recombination on the re- 
maining parameters is shown in Figure [S] where we 
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FIG. 6: Likelihood distribution function from WMAP5 on 
the optical depth (top panel), the baryon density (center 
panel) and the cold dark matter density (bottom panel). 
The red dashed line are the results assuming standard re- 
combination, the solid black line are under the hypothesis 
of delayed recombination. 



plot the 1-D likelihood functions for the optical depth 
T and the baryon and cold dark matter densities Wb 
and iOc derived from the WMAP5 data with and with- 
out standard recombination. While the changes in 
the optical depth are minimal since this parameter 
is mainly fixed by large scale polarization data, de- 
layed recombination may slightly bias the current con- 
straints towards smaller wj, and larger ujc respect to 
the standard case. 

Another parameter worthy of study, and on which 
the CMB usually provides clear and uncorrelated con- 
straints, is the shift parameter, TZ, (see eg. [4l[) which 
defines the angular scale of the acoustic oscillations. 
The shift parameter is generally considered a simple, 
geometric way to implement information from CMB 




FIG. 7: Constraints on the z^-ea plane from WMAP5 
(empty contours) and WMAP5+ACBAR (filled contours). 
A degeneracy along the relation Zt{ea) = 2;*(0)(1 -I- 
3e<,)"°°''^ is evident. 



in determining constraints from independent data sets 
such type la supernovae or galaxy clustering. How- 
ever this parameter is not directly measured by CMB 
observations but is essentially a byproduct of the cos- 
mological parameter inference based on CMB spec- 
tral data analysis, and most crucially the recombina- 
tion history. Therefore, using the shift TZ to constrain 
cosmological parameters in combination with other 
datasets may suffer from model dependent assump- 
tions. It is subsequently important to investigate the 
effects modified recombination can have on this key 
parameter. We find WMAP5 with delayed recombi- 
nation constrains TZ = 1.737 ±0.028 at 68% cd., to be 
compared with 7^ = 1.713± 0.020 at 68% c.l. for stan- 
dard recombination. Such a 1 — cr shift in 7?. should 
not, however, alter in a significant way the current 
constraints on dark energy parameters derived from 
combined analyses assuming the standard recombina- 
tion. 

As expected, modified recombination strongly af- 
fects the constraints on the redshift of recombination, 
z*. According to the latest WMAP5 results, the red- 
shift of recombination z* is z* = 1090.51 ±0.95 with a 
better than 0.1% precision. However this impressive 
result is mainly driven by the assumption of standard 
recombination itself. With delayed recombination, we 
found z* = 1078.2 ± 10.9 at 68% c.l. from WMAP5 
alone and z* = 1084.2 ± 9.9 with WMAP5+ACBAR, 
so that the uncertainty in the measurement of z* is 
increased by an order of magnitude. 

In Figure [7] we plot the 2-D likelihood contours 
for z* and e^. A degeneracy following the relation 
z*(e„) ~ z*(0)(l + 3ec<)-°°''2 (see [13) is clearly evi- 
dent and explains why the constraints on z* are much 
broader in the presence of delayed recombination. 
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FIG. 8: Constraints on the Q^j-Qa plane in the presence 
of delayed recombination from WMAP5+ACBAR (empty 
contours) and after adding an external prior Qj^j = 0.30 ± 
0.05 (filled contours). 



B. Modified recombination with curvature 

As mentioned in the previous section, delayed re- 
combination decreases z^,, increasing the size of the 
sound horizon at recombination. This shifts the po- 
sition of the acoustic peaks in the CMB spectrum 
as ^ ~ (^^/^jtandard)i/2^ ^g therefore expect a de- 
generacy between Eq. and the curvature ujk: namely a 
lower recombination redshift shifts the CMB peaks to 
larger angular scales (smaller ^'s) yieldin g op en mod- 
els more consistent with the data (see [13] )■ How- 
ever, as already pointed out in [l3|, measurements of 
the multiple peaks in the angular spectrum should 
break this degeneracy. In this section we perform an 
analysis of delayed recombination in curved models 
and check if a flat universe is still consistent with the 
data and test to what extent open models can be in 
agreement with observations. Including curvature, we 
found that the constraints on WMAP5+ACBAR are 
slightly relaxed with e^ < 0.39 and e^ < 0.058 at 
95% c.l. However the effect on curvature from de- 
layed recombination is small. Marginalizing over the 
recombination parameters we get Uk = — O.OSS^q^oo 
(WMAP5-hACBAR) at 95% c.l.. to be compared with 
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obtained using the 



same datascts and priors but with standard recombi- 
nation. The shifts towards open models due to delayed 
recombination is therefore at the level of few percent. 
When external priors are added, we see little vari- 
ation in the current constraints on Q\ in a curved 
universe. We show this in Figure [S] where we plot 
the usual constraints on the Qm-^a but including de- 
layed recombination. The changes are minimal re- 
spect to standard recombination. Adding a prior on 
the matter parameter J7m = 0.30 ± 0.05 yields a 
constraint J^a = 0.70 ± 0.04 (to be compared with 
r^A = 0.712 ±0.037 with standard recombination) still 



suggesting the presence of a cosmological constant at 
high significance. 



C. Modified recombination and massive 
neutrinos 

Cosmological neutrinos have a relevant impact on 
cosmology since they change the expansion history of 
the universe and affect the growth of perturbations. 
CMB anisotropies can constrain neutrino masses in- 
directly since variations in the gravitational potential 
change the shape of the small scale CMB angular spec- 
trum. 

The recent WMAP5 data have provided an upper 
limit to the sum of neutrino masses of Ettt,^ < l.SeV 
at 95% c.l. (H). Since this is the best upper limit 
current available on absolute neutrino masses it is 
certainly timely to investigate how this bond could 
change when delayed recombination is included. We 
have found that from WMAP5-f ACBAR the bound 
on neutrino masses is Tirrii, < 1.26ey at 95% c.l. un- 
der delayed recombination larger than Etti^ < 1.13eT^ 
obtained under standard recombination with the same 
dataset and priors. We can therefore conclude that 
delayed recombination increases the current bound by 
~ 10%. 



D. Forecasts for the Planck Surveyor 

In this section, we forecast the future constraints 
achievable on modified recombination from the Planck 
satellite experiment. We have constrained cosmic 
parameters assuming simulated Planck mock data 
with a fiducial model given by the best fit WMAP5 
model (with standard recombination) and experimen- 
tal noise described by 
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with w ^' ^ = 63^K as the temperature noise level (we 
consider a factor \'2 larger for polarization noise) and 
^FWHM = 7' for the beam size. We take fsky = 0.65 
as sky coverage. 

With this configuration, we found the following con- 
straints from Planck: e^ < 0.01 and Ei < 0.0005 at 
95% C.L, representing more than an order of magni- 
tude in improvement in ea and two order of magnitude 
improvement in e^, respectively, in comparison to cur- 
rent WMAP5+ACBAR constraints. 

In Figures [S] and [TO] we show the effects of delayed 
recombination on constraining cosmological parame- 
ters from Planck. As we can sec, the precision small 
scale temperature measurements of Planck markedly 
reduce the Ug — t — Ci/ca degeneracies, e^ is further 
constrained by small scale polarization measurements. 
Subsequently a modified recombination scheme should 
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FIG. 9: Likelihood distribution function from a Planck-like 
satellite experiment (see text for the experimental configu- 
ration) on the scalar spectral index (top panel), the baryon 
density (center panel) and the cold dark matter density 
(bottom panel). The red dashed line are the results as- 
suming standard recombination, the solid black line are 
under the hypothesis of delayed recombination. 



result in values of the spectral index Ua being only 
slightly skewed towards more positive values in com- 
parison to the standard scheme, and the Planck satel- 
lite will therefore be clearly able to discriminate n^ 
from a HZ spectrum with high significance even in 
the case of delayed recombination. The significantly 
larger uncertainties in ip and Hq currently introduced 
if delayed, rather than standard, recombination is con- 
sidered with the WMAP5 -hACBAR data would be 
significantly curtailed with Planck measurements. 
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FIG. 10: 2D likelihood contour plots at 68% and 95% c.l. 
from a Planck-like satellite experiment (see text for the 
experimental configuration) for the Hubble parameter (top 
Panel) and age (Bottom Panel) versus the rcombination 
parameter £„. 



IV. CONCLUSIONS 

In this paper, we have updated the upper bounds 
that can be placed on the contribution of extra Ly- 
a and ionizing photon-producing sources in light of 
the new WMAP and ACBAR data. We have found 
that, adopting a simple parametrization using con- 
stant, effective values for e^ and e^, the WMAP5 data 
constraints £„ < 0.31 and e^ < 0.053 at the 95% level. 

We have studied the implications of delayed recom- 
bination on several current cosmological constraints 
derived from CMB anisotropics. While we found no 
relevant changes on current CMB estimates of cur- 
vature, optical depth, baryon and cold dark matter 
densities, delayed recombination proves to have a sig- 
nificant impact on current constraints for infiationary 
parameters. Current conclusions for theoretical infla- 
tionary models, motivated by evidence for a deviation 
from scale invariance from recent CMB data in the 
standard recombination scheme, would have to be re- 
laxed when delayed recombination is considered, since 
the inflationary spectral index of scalar perturbations 
is in complete agreement with a scale invariant HZ 
spectrum. 

The delayed recombination also has an impact on 
estimates of the shift factor, 7?., and recombination 
redshift, z*, currently used in estimating angular di- 



ameter distance and baryon acoustic oscillation mea- 
sures of the cosmic expansion history. Moreover, con- 
straints on particle physics parameters like the neu- 
trino mass can also be relaxed when non-standard re- 
combination is considered. 

Physically motivated models for non-standard re- 
combination, like those based on primordial black hole 
or super heavy dark matter decay, arc possible and 
provide a good fit to the current data. Future ob- 
servations in both temperature and polarization, as 
those expected from the Planck satellite will be needed 



if we are to more stringently test and reduce the de- 
pendency of other cosmological parameters on delayed 
recombination models. 
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